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Abstract 

It is pointed out that electrostatic and electromagnetic waves prop¬ 
agating in perpendicular direction to the density gradient Vno = 
x|^r| and external magnetic held Bq = Bqz in hybrid frequency range 
Qi -C OJ <C (where ilj = and j = e, i ) have several applications 
in space and laboratory plasmas. The electron plasma waves are im¬ 
portant for the future experiments on heavier ion plasmas as well as 
the results are applicable to the experiments performed to create pure 
pair ion fullerene plasmas. The electrostatic modes are shown to exist 
in upper F-region of terrestrial ionosphere and electromagnetic waves 
are applied to cylindrically bounded heavier ion plasmas. 


1 Introduction 

The electrostatic drift type waves in hybrid frequency range (H* < oj < fi e ), 
where Qj = is the gyro frequency of j-th species and j = (e, i ) have been 
investigated in the recent past (Eliasson et al. 2006; Saleem and Eliasson 
2011) in the light of future plans to perform experiments with heavier ions. 
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The study of heavier ion plasmas is also important due to their relevance 
to experiments performed to produce pair ion fullerene plasmas (Oohara 
and Hatakeyama 2003; Oohara et al. 2005; Oohara eta al. 2009) which 
have generated controversy among theoreticians (Saleem 2006; Verheest 2006; 
Kono eta al. 2014). 

In heavier ion plasmas the gap between cyclotron frequencies of ions O,; 
and electrons Vt e is large and hence some waves may appear in the frequency 
range 0* <C uj < fl e and for these waves the ion dynamics can be ignored. 
However, for the space plasmas of hydrogen or fusion plasmas of deuterium 
(D)-tritium (T), the role of ions is important in the hybrid frequency waves 
due to small ratio of masses —. 

rtii 

More than a decade ago (Huba 1991; Huba eta al. 1995), a purely trans¬ 
verse wave, the magnetic drift wave (MDW), was proposed to exist in inho¬ 
mogeneous density electron plasmas with linear dispersion relation 

w = kyVAi ^-) = A lk*(£n e ) (1) 

\LilJn <Vy 

where A e = — is the electron skin depth. u ve = / 4nn o e ~ j s the elec- 

(Ope. * y ITle 

tron plasma frequency, L n = 0- is the density gradient scale length, and 
va = (i s Alfven speed. The density gradient is assumed to be in 
the positive x-direction Vno = x \ \, B 0 — zB 0 and wave propagates 

in y-direction perpendicular to both density gradient and external magnetic 
field. The linear dispersion relation (1) was derived using electron magne¬ 
tohydrodynamics (EMHD) equations. In standard EMHD model (Kingssep 
et al. 1990), the electrons are considered to be inertia-less (m e —> 0) and 
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ion dynamics is ignored in the limit (m, —> oo) for time scale r such that 
fl” 1 , Wpg 1 Cr< 12” 1 , Wpi- 1 that is for frequency range 12;, ca p ; CwC uj pe , 12 e . 

Recently (Saleem 2014), it has been proposed that one should use equa¬ 
tion of motion to determine the perpendicular electron velocity components 
v e = ( v ex ,v ey ) instead of using Ampere’s law for v e . If electron velocity is 
determined appropriately using equations of motion, then one notices that 
compressibility cannot be ignored in the MDW model. The wave becomes 
partially transverse and partially longitudinal and in the presence of density 
gradient its dispersion relation is modified as 



In heavier ion plasmas, this wave can play an important role and it has 
not been investigated in such plasmas so far. In lighter element plasmas, the 
frequency window 12; < u < 12 e is narrower and this wave couples with the 
lower hybrid oscillations (Saleem 2014). The electrostatic (ES) lower-hybrid 
wave (LHW) is important due to its role in laboratory plasma heating and 
spacecraft observations of terrestrial ionosphere and magnetosphere . In 
the presence of density gradient, the LHW couples with the electrostatic 
electron convection mode (ECM) (Yu et al. 1985). Both the LHW and 
ECM propagate obliquely making a small angle with the external constant 
magnetic field satisfying k z k± and 12; < u < 12 e while k z (k±) are parallel 
(perpendicular) components of the wave vector, and u is the wave frequency. 

The pure transverse electromagnetic wave known as the ordinary O- mode 
(Chen 1984) in plasmas couples with the transverse current and its linear 
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dispersion relation becomes ( u 2 = u 2 e + c 2 k 2 ) . The extraordinary wave 
(X-rnode) is also an electron time scale wave and it is partially transverse 
and partially longitudinal (Chen 1984). Both 0 and X modes are high fre¬ 
quency and short wavelength waves and do not need ion dynamics for their 
description, in general. 

In this work, it is pointed out that the extraordinary (X-) mode in the low 
frequency limit becomes drift type wave in an inhomogeneous density plasma 
which can couple with the lower hybrid oscillations when the ion dynamics is 
taken into account. This wave propagates in the direction perpendicular to 
the external magnetic field and has frequency in the lower hybrid frequency 
range fl* <C u <C kl e . The basic wave is described in cartesian geometry and 
its electric field structure is found in a cylindrically bounded plasma by solv¬ 
ing the coupled Poisson equation and Ampere’s law. In heavier ion plasmas 
the ion dynamics may be ignored. 

Here we mainly discuss the electromagnetic extraordinary drift wave in cylin¬ 
drically bounded plasma for its application to future experiments on heavier 
ion and pair ion fullerene plasmas. Then the eigenvalue problem is also solved 
in a cylindrically bounded plasma with Gaussian density distribution. 


2 Extraordinary Drift Wave 


Let us consider the extra-ordinary wave with Bq 
E = xE x + yE y in an inhomogeneous plasma Vn 0 = 


= B 0 z, k = k y y and 

„, dn 0 


—x\ 


dx 


The Poisson 
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equation yields 


O 2 


~Upe~ 


K n o 1 Z7- ; L0 P e (O 

]T-^. - (Si 


e + "7 “ w )^1j: — 0 (3) 

ky 


and Maxwell’s equation provides us with another relation between E\ x and 
Ely as follows 


l(w 2 - oV)(l - % 


oo 


^pe 


a 


CO 


+ uJL( 


^ + < 4 > 


If we neglect K n terms, then (3) and (4) are similar to the set of equations 
(4-101) of (Chen 1984) which gives dispersion relation for X-mode as 


/ 2 2\[ 2 2 2; 2/1 ^e\l t^V ^ e \2 

{U - U H )[U - U H - C k (1 - —)] = (- ' 


CO 


00 


( 5 ) 


where co 2 H = (c o 2 e + fig) is the upper hybrid frequency wave. The X-mode 
is well studied in literature but the low frequency regime of this wave in 
the presence of density gradient has not been investigated. In lighter ion 
plasmas, the wave frequency becomes closer to 0,: or {kl e Vti) 1 ^ 2 in the limit 
oo 2 <C oo 2 e , fig , therefore the role of ion dynamics must be included. In this 
case Eq. (2) becomes, 
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Long ago (Yu et al. 1985), the electrostatic electron convection mode was 
investigated which has the dispersion relation 


where 


rru A ; 2 , 

" = tJ ” 11 + ^if 1 


UJ V ( JxyL n 


( 7 ) 


( 8 ) 


for k z <C k y . 

The coupled dispersion relation of electron convection mode and lower hybrid 
wave was presented in equation (5) of Ref. |13] as, 


k 


, .2 2 

| _ yV , 

w 2 k 2 [ L n ojQ e 


k 2 E 1 

— H —-1 = 0 
2 ^ fiaJ 


a; 


( 9 ) 


Yu et al. (1985) discuss the obliquely propagating two modes only; the lower 
hybrid wave and electron convection mode. 

It is pointed out here that the Eq. (9) also contains an electron mode which 
propagates purely in the perpendicular direction to the external magnetic 
held and density gradient. 

Let k z = 0, k x = 0, and put uj pi = 0 (i.e. ions stationary), then Eq.(9) 


reduces to 


1 


for 


UJ = 


Q 2 s v k. 


n e ) 


( 10 ) 


(i + -f) 


e 

a; 2 

pe 


In heavier ion plasmas (like Barium or Fullerene), the above wave will be 
relevant. In lighter element plasmas, this will couple with the lower hybrid 
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frequency and Eq.(lO) will become, 


(1 H-^")^ 2 + T ~^e — ~ 0 (11) 

W pe ky 

Both dispersion relations (10) and (11) have been presented by Saleem (2014) 
but it seems important to point out that these simple dispersion relations 
were present in earlier literature implicitly. 

3 Cylindrically Bounded Plasmas 

The waves discussed in previous section have different regimes of application. 
The electrostatic modes have good applications both in space and laboratory 
plasmas, while the electromagnetic modes are more relevant to plasmas with 
1 < A Iky and this limit is generally satisfied in laboratory plasmas. 

The heavier ion fullerene plasmas have been produced in a cylinder of length 
l = 90 cm and diameter Rq = 3 cm (Oohara et al 2005). The planned Barium 
plasma experiments will probably be performed in a cylindrical geometry (Q 
Machine). 

Therefore, in this section, we present the form of local modes in cylindrical 
geometry and plasma is assumed to have Gaussian density profile, 

no(r) = N 0 e~ r2/R o (12) 

where i?. 0 is the cylinder radius and N 0 is the density at the axis. The external 
magnetic field is along the 0 -axis and waves propagate in d-direction. The 
external magnetic field is along 0 -axis, the perturbed electric and magnetic 
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fields are, respectively, E = rE r + OEg and B = zB. The perturbations are 

assumed to have the form i/j = •0(r)e t ( m6l_a,t ) where m is the azimuthal wave 

number and kg — —. 

r 

Poisson equation takes the form, 


(Li - ^ )E, 


m 


— L - 

m 


)d r Eg 


T 2 

{-uQ e + ^ + {to 2 - nM—LE r + Ll—ld r E r } (13) 
m m m 


where 


4 = 

Q 2 

Lq = (1 H- — ■ 

(jJ L 

pe 

In the limit u <C f2 e , cUp e , Ampere’s law gives, 

E r = 2 K (T” + ~i rn ) E o + <9 r ^} (14) 

A 2 m 2 \l e r 1 r 


a Local modes in cylinder 

The density profile (12) gives in cylindrical geometry K n 
the local wave dispersion relation takes the form, 


2 r 

- —, therefore 
Rq 


{1 + (1 + 


O 2 X 2 

* L e \ A e 2 a 2 , /o 

—)—m \(jj + (2m 


y 
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)f2 e o; — 0 



= 0 


(15) 



which in a heavier ion plasma with static ions becomes 


c o = 


(2 m^)Q e 

K 0 


02 \ 2 

{ 1+(1 + ^)^ 2 } 

uj z r z 

pe 


(16) 


In electrostatic limit, Eq. (15) reduces to 


Q 2 r 2 

(1 + -§-)(, 0 2 + 2 = 0 


cu 


pe 




(17) 


which in a heavier ion plasma gives, 


u = 


(2m 2 )^e 

o 2 

(l + -f) 


cu. 


pe 


(18) 


b Eigen Modes of X - drift wave 

Here we analyze the low frequency W-mode in a cylinder with small variation 
of amplitudes along radial direction assuming that the fields vanish at the 
boundary r = Rq. We use the conditions 


m 


K n d r E e < E d , \K n d r E e \ < |— d r E e \ <m \E t 


(19) 


and consider the waves in the region 0 <C r < Rq. 

Then the Poisson equation and Ampere’s law after simplification yield, 

r 2 d 2 E e + rd r Ee + (o- 2 r 2 — m 2 )E g = 0 (20) 
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where 


a = 


, u „ m ^ 

— + 2^n e 


Q 2 

— (i h— )u 2 } 


w, 


pe 


The r-dependent general solution of Eq. (20) is 


E e (r) = .4, J m (ar) + A 2 Y m (ar ) 


( 21 ) 


( 22 ) 


where J m , are Bessel functions of order m and are constant am¬ 

plitudes. 

At r = 0, Y m blows up, therefore we put = 0 and the physical solution 
becomes, 

E e (r) = A\J m {ar) (23) 

We choose the boundary condition Eg — 0 at r — R 0 , therefore J m (ai? 0 ) = 0 
gives k a = aRo where k a are the zeros of the Bessel function J m . 


4 Applications 

First the electrostatic waves are shown to exist in F-region of terrestrial 
ionosphere considering Uq = 7.5 x 10 3 cm -3 , B 0 = 0.3 Gauss in hydrogen 
plasma assuming the density gradient scale length L n = 5 x 10 4 cm = 500 
m (i.e., K n — 2 x 10 5 cm _1 ) . This plasma has Vt e = 5.2 x 10 6 rads _1 , fl; = 
2.8 x 10 3 rads _1 , A e = 6.1 x 10 3 cm, ojlh = 1-2 x 10 5 rads _1 and we choose 
k y : (0.005 — 0.05) cm so that 1 <C k 2 y \ 2 e holds and hence only electrostatic- 
waves are investigated within local approximation K n <C k y . In applications 
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we shall use Vno = x | ^ | i.e. density gradient along + ve x-axis and then u 
of Eqs. (10) and (11) will be positive. 

The Fig. (la) shows u vs k y corresponding to electron plasma dipersion 
relation Eq. (10). It is obvious from Fig. (la) that hh < u <C kl e and hence 
ion dynamics should be considered. 

In Fig. (lb), the plot of u vs k y of Eq. (11) is shown and flj < cu < Vt e 
is satished. The wavelength of the wave X u — : (12 — 1) m satisfies the 

condition A w <C L n . For small k y : (10 -3 — 10 -2 cm -1 ) , the wavelength 
becomes larger and Fig. (2a) and Fig. (2b) show, respectively, the plots u 
vs k y of Eqs. (10) and (11). Corresponding to k y = 10 -3 , the wavelength 
becomes A w = 60 m and still A w <C L n holds. 

The dispersion relations for electromagnetic waves in pure electron and electron- 
ion hydrogen plasmas are shown in Fig. (3a) and Fig. (3b), respectively us¬ 
ing Cartesian geometry and JET parameters i.e. n 0 - 2 x 10 14 cm -3 , B 0 = 
10 4 Gauss, and K n = 0.1cm -1 (L n = 10 cm). In this case, = 1.7 x 
10 n rads -1 , Vt e = 1.7xlO n rads -1 , oj lh = 4.1 x 10 9 rads -1 and A e = 0.0375 cm. 
Since u is closer to fi*, therefore Fig. (3a) does not satisfy the static ions 
condition -C u. The ion dynamics is needed and hence dispersion relation 
of Eq.(6) is shown in Fig. (3b). 

In case of JET parameters, the cylindrical coordinates are better to use. 
However, in this case, we take i?o as average value of JET minor radius 
R 0 = 1.67 m. The Fig. (4a) shows the plot cavsr of Eq. (16) corresponding 
to m = 500 (dashed curve) and m = 1000 (solid curve). This EM wave is 
relevant for 0 <C r away from axis. Near axis, u is closer to lower hybrid 
frequency, therefore ions should be considered to be dynamic as shown in 
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Fig. (4b). Here we obtain uj < fl*. Therefore, the dispersion relation of Eq. 
(16) does not seem to be valid for H-plasma in this case. 

Now we consider Barium plasma in a cylinder with n 0 — 2 x 10 13 cm -3 and 
B 0 = 10 4 Gauss for L n = 10 cm. In this case, = 6.9 x 10 5 rads -1 , and we 
choose R() = 100 cm. 

The Fig. (5a) shows u vs r for Eq. (16) and heavier ions are considered to 
be stationary because Q,; <C w holds. However, the dispersion relation is not 
valid near r = 0, the singular point. 

In Fig. (5b), u vs r from Eq. (15) is plotted for Barium plasma taking into 
account the ion dynamics. 

Finally, in Fig. (6), we present the plot of Eq. (23) the solution of eigenvalue 
Eq. (20) for Eg in a cylinder of Rq = 100 cm for JET parameters just as 
an illustration taking A\ = 1 and choosing u = 3 x 10 9 s _1 from Fig. (3b). 
This figure shows how the amplitude of rapidly oscillating Eg with frequency 
w = 3x 10 9 s _1 decays in the cylindrical (Barium) plasma for 40 < r < Ro 
corresponding to m = 500. 

5 Discussion 

The applications of electrostatic and electromagnetic high frequency drift 
waves have been discussed with reference to upper F-region of terrestrial 
ionosphere and cylindrically bounded laboratory plasmas. These waves prop¬ 
agate in the direction perpendicular to the external magnetic held and the 
density gradient. The waves fall in frequency range of lower hybrid waves 
(LHW). The LHW is electrostatic and propagate obliquely making a small 
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angle with the ambient magnetic held with fey <C k±, but the waves discussed 
in this work have k\\ = 0. 

A great deal of literature exists on waves in inhomogeneous density plas¬ 
mas and most of the work has been done on the low frequency drift waves 
(u <C O,:) both in electrostatic and electromagnetic limits due to their ap¬ 
plications in particle and energy transport. A few authors have reported 
two new modes in the LH frequency range; one is the electrostatic electron 
convection mode (Yu et al. 1985) and the other is a purely transverse wave; 
the magnetic drift wave (MDW) (Huba 1991; Huba et al. 1995) which has 
k\\ = 0. In a recent paper (Saleem 2014), it has been presented that the com¬ 
pressible effects should be added to MDW and then the dispersion relation 
will be modified and one finds two different limits of such electromagnetic 
drift waves; one is applicable to pure electron plasma and the other contains 
the effect of lower hybrid oscillations. 

Here we have pointed out that the electromagnetic waves discussed in by 
Saleem (2014) are actually the limiting case of the well-known high frequency 
extraordinary (X-) wave. Since X-mode is believed to be a high frequency 
short wavelength mode, probably that is the reason that density gradient 
effects on this mode are considered to be small. But it is interesting to look 
at a low frequency limit of this mode which gives two important dispersion 
relations; one is applicable to purely electron plasmas (Eq.(2) )and the other 
contains the ion dynamics (Eq. (6)) as well. 

The electrostatic waves have been shown to exist in upper F-region iono¬ 
sphere where 1 <C A 2 e k^ holds. The elctromagnetic waves have been applied 
to a cylindrically bounded plasma with Gaussian density distribution. For 
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illustration point of view, the JET plasma parameters have been used for 
which A Iky < 1 is satisfied. 

The waves discussed here with their applications to both space and labo¬ 
ratory plasmas can have many more applications. The further investigation 
of these waves will be useful for future laboratory experiments particularly 
with heavier ion plasmas like barium and fullerene. They will also be appli¬ 
cable to space plasma observations. 
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List of Figures 


1. Fig.(la) The plot to vs k y of Eq. (10) shows that to is closer to £1* and 
hence in hydrogen like plasmas, the ion dynamics is necessary to be 
considered. 

2. Fig.(lb) The to vs k y is plotted using Eq. (11) for dynamic ions of 
H-plasma. 

3. Fig.(2a) Plot of Eq. (10) for larger wavelengths still keeping <C L n . 

4. Fig.(2b) Plot of Eq. (11) for larger wavelengths still keeping <C L n . 

5. Fig. (3a) The plot for to vs k y with JET parameters show fl* < to and 

therefore in H-plasma Eq. (2) is not valid. 

6. Fig. (3b) For H-plasma with JET parameters to vs k y is plotted with 
dynamic ions (Eq 6) using Cartesian geometry. 

7. Fig. (4a) The plot of pure electron plasma for H-case Eq. (18) as to vs 
r for m = 500 (dotted line) and m = 2000 (solid curve) corresponding 
to JET parameters. Obviously to is closer to ! !,. 

8. Fig. (4b) Here to v r is plotted considering H-ions to be dynamic corre¬ 
sponding to JET parameters and we find to < fli . Therefore, Eq.(6) is 
not valid for H-plasma. 

9. Fig. (5a) Eq. (2) is plotted for to vs r for Ba plasma and in this case 
ions can be considered to be stationary because Oj holds. 
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10. Fig.(5b) The electromagnetic dispersion relation of Eq. (6) is plotted as 
u> vs r for Ba plasma considering ions to be dynamic and fl* < u < fl e 
holds. 

11. Fig. (6) The radial dependence of rapidly oscillating Eg of Eq. (23) 
corresponding to u = 3 x 10 9 s -1 of Fig. (3b) in cylindrically bounded 
inhomogeneous plasma. 
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Usman et. al. Fig. (8£) 






























